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A particle of mass /i moves on a circular orbit of a nonrotating black hole of mass M. Under 
the restrictions ji/M <C 1 and v <^ 1, where v is the orbital velocity (in units in which c = 1), we 
consider the gravitational waves emitted by such a binary system. The framework is that of black- 
hole perturbation theory. We calculate E, the rate at which the gravitational waves remove energy 
from the system. The total energy loss is given by -B = -E°° + , where E°° denotes that part 
of the gravitational-wave energy which is carried off to infinity, while E^ denotes the part which is 
absorbed by the black hole. We show that the black-hole absorption is a small effect: E^ /E ~ n*. 
This is explained by the presence of a potential barrier in the vicinity of the black hole: Most of 
the waves propagating initially toward the black hole are reflected off the barrier; the black hole is 
therefore unable to absorb much. The black-hole absorption (and indeed any other effect resulting 
from imposing ingoing- wave boundary conditions at the event horizon) are sufficiently small to be 
irrelevant to the construction of matched filters for gravitational-wave measurements. To derive 
this result we extend the techniques previously developed by Poisson and Sasaki for integrating the 
Regge- Wheeler equation. The extension consists of an explicit consideration of the horizon boundary 
conditions, which were largely ignored in the previous work. Finally, we compare the wave generation 
formalism which derives from perturbation theory to the post-Newtonian formalism of Blanchet and 
Damour. Among other things we consider the corrections to the asymptotic gravitational- wave field 
which are due to wave-propagation (tail) efi'ects. The results obtained using perturbation theory 
are identical to that of post-Newtonian theory. 
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I. INTRODUCTION AND SUMMARY 

A. Gravitational waves from coalescing compact 
binaries 

Coalescing compact binary systems, composed of neu- 
tron stars and/or black holes, have been identified as 
the most promising source of gravitational waves for 
kilometer-size interferometric detectors 

The construction of three such detectors should be 
completed by the turn of the century. The American 
LIGO (Laser Interferometer Gravitational-wave Obser- 
vatory) project jsj involves two 4 km detectors, one sit- 
uated in Hanford, Washington, the other in Livingston, 
Louisiana; construction has begun at the Hanford cite. 
The French-Italian Virgo project jj] involves a single 3 
km detector, to be built near Pisa, Italy. 

According to recent estimates , these interferome- 
ters should be able to detect approximately three binary 
coalescences per year, at distances of order 200 Mpc. If 
the detectors achieve the so-called "advanced detector" 
sensitivity then the limiting distance will be of the 
order of 1 Gpc, and the event rate will be increased by a 
factor of order 100. 

The emission of gravitational waves by a compact bi- 
nary system causes the orbits to shrink and the orbital 
frequency to increase The gravitational- wave fre- 
quency is given by twice the orbital frequency; it also 
increases as the system evolves. The LIGO and Virgo 
detectors are designed to operate in the frequency band 
between approximately 10 Hz and 1000 Hz g. This cor- 
responds to the last several minutes of inspiral, and pos- 
sibly the final coalescence (depending on the size of the 
masses involved) , of a compact binary system Q . 

During this late stage of orbital evolution the orbital 
velocity is large {v/c ranging from approximately 0.07 to 
0.35 if the binary system is that of two neutron stars [^), 
and so is the gravitational field [the dimensionless gravi- 
tational potential is of order {v/cf']. Accurate modeling 
of the inspiral must therefore take general-relativistic ef- 
fects carefully into account. On the other hand, other 
complications, such as tidal interactions (important only 
during the last few orbital cycles ||^,|l^) and orbital ec- 
centricity (reduced to very small values by gravitational 
radiation reaction 12|) can be safely ignored. 

The gravitational waves emitted by a coalescing com- 
pact binary carry information about the source — the 
wave forms depend on the parameters describing the 
source. These include the orientation and position of 
the source in the sky, its distance, and the masses and 
spins of the companions. A major goal in detecting these 
waves is to extract this information [p][2[jT^. 

Because the wave forms can be expressed, at least in 
principle, as known functions of the source parameters, 
the technique of matched filtering can be used to esti- 
mate the values of these parameters [p]j7p^ : One cross- 
correlates the gravitational-wave signal with a set of theo- 



retical wave forms (templates); the template with param- 
eters closest to the actual values maximizes the signal-to- 
noise ratio. For this technique to work, the gravitational- 
wave signal must be accurately predicted in the full range 
of frequencies described above. It can be estimated that 
the waves oscillate a number of times of order IC* as the 
frequency sweeps from 10 Hz to 1000 Hz Q. To yield 
a large signal-to-noise ratio the template must stay in 
phase with the measured signal. The accurate prediction 
of the phasing of the wave is therefore a central goal for 
theorists [7[p^; an accuracy better than one part in 10* 
is required. 

Gravitational radiation reaction causes /, the 
gravitational-wave frequency, to increase as the system 
evolves. The phasing of the wave is determined by /(/), 
the rate of change of frequency as a function of frequency. 
If the orbital energy E is known as a function of /, then 
E{f) determines the phasing. Below we will assume that 
the relation E{f) is known (or can be calculated), and 
focus on the energy loss. 



B. Energy loss: post-Newtonian theory 

The energy radiated by a compact binary system can 
be calculated, approximately, using post-Newtonian the- 
ory ||]. The post-Newtonian approximation is based 
upon an assumption about the smallness of the orbital 
velocity. In units in which c = 1 (we also put G = 1), we 
demand f ^ 1. On the other hand, nothing is assumed 
about the mass ratio: If M denotes total mass and /i 
reduced mass, then /x/Af is allowed to cover the whole 
range (0, 1/4]. For the problem under consideration v is 
small only during the early portion of the inspiral, when 
/ is near 10 Hz. It is therefore expected that calculations 
will have to be pushed to extremely high order in v. An 
intensive effort to do just this is now underway |8j. 

The phrase "post-Newtonian theory" is used here 
loosely. A wave-generation formalism suitable for 
high-order calculations was developed by Blanchet and 
Damour [p"5|-^. This combines a post-Newtonian ex- 
pansion (based, in effect, on the assumption c cxj), for 
the region of spacetime corresponding to the near zone, 
with a post-Minkowskian expansion (based, in effect, on 
the assumption G — > 0), for the region of spacetime out- 
side the source. The external metric is matched to the 
near-zone metric, and the radiative part of the gravita- 
tional field is extracted. The radiative field is charac- 
terized by an infinite set of multipole moments; the re- 
lation between these radiative moments and the source 
moments is revealed by matching. 

The rate at which energy is radiated can be calcu- 
lated from the radiative moments, and expressed as an 
expansion in powers of v. (In fact, the expansion also 
involves Inu at high order jljjl^,|l^.) Thus far, the en- 
ergy loss has been calculated accurately to order t>*, or 
post^-Newtonian order. Schematically, 
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where Eq^f = {?,2/b){^i/MYv'^° is the leading-order 
(Newtonian) quadrupole-formula expression Here 
and below, we define v = {MVlYl'^ = (ttM f)^'^, where 



is the a ngul ar velocity. 



In Eq. (1.1) the 0(v ) term was first calculated by 



Wagoner and Will |^0| and is due to post-Newtonian 
corrections to the equations of motion and wave gener- 
ation. The 0{v^) term comes from two different effects: 
(i) tail terms in the waves (the gravitational waves scat- 
ter off the spacetime curvature as they propagate from 
the near zone to the far zone), first calculated by Poisson 
pl| using perturbation theory (see below), and subse- 
quently by Wiseman and Blanchet and Shafer p^ ] 
using post-Newtonian theory; and (ii) spin-orbit interac- 
tions (present if at least one of the companions is rotat- 
ing), first calculated by Kidder, Will, and Wiseman [ p4| . 
The 0(y'^) term originates from spin-spin interactions (if 
both stars are rota ting ), as was calculated by Kidder, 
Will and Wiseman [Q, and from post^-Newtonian cor- 
rections to the equations of motion and wave generation, 
as was recently calculated by Blanchet, Damour, Iyer, 
Will, and Wiseman |2^. An additional contribution to 
the 0(w"') term will be present if at least one of the stars 
has a nonvanishing quadrupole moment; this effect has 
not yet been calculated, except for the special case con- 
sidered by Shibata et al. 126| (see below). 



tensor T"'^ which creates a perturbation in the gravita- 
tional field. If /z <C M the Einstein field equations can 
be solved perturbatively about the black-hole solution. 
The perturbations propagate as gravitational waves. By 
solving the perturbation equations one can calculate the 
wave forms, as well as the energy and angular momen- 
tum transported. The energy and angular momentum 
are partly carried off to infinity, and partly absorbed by 
the black hole. 

The motion of the particle must be specified before the 
perturbation equations are integrated. (The particle's 
stress-energy tensor, which acts as a source term in the 
equations, is a functional of the world line.) In general, 
the motion is affected by the loss of energy and angular 
momentum to gravitational waves — there is a radiation 
reaction force |pQ]. However, the radiation reaction can 
be neglected to lowest order in /i/Af , and the motion can 
be taken to be a geodesic of the black-hole spacetime. 
The stress-energy tensor is then completely specified. 

Because the source term in the perturbation equations 
is known exactly, to integrate these equations amounts 
to solving an equation for wave propagation in curved 
spacetime. (In post-Newtonian theory, wave propaga- 
tion is only part of the problem. One must also solve 
for the near-zone physics.) For large velocities the wave 
equation must be integrated numerically. When com- 
bined with a slow-motion approximation, the equations 
of perturbation theory can be integrated analytically. 



C. Energy loss: perturbation theory 



D. Perturbation theory: a survey 



Post-Newtonian theory is based upon a slow-motion, 
weak-field approximation. Because the orbital velocity is 
large during the late stages of the inspiral, and because 
the post-Newtonian expansion is not expected to con- 



verge rapidly {if at all; the coefficients in Eq. (1.1) are 



typically large and grow with the order [|18|}, it is not 
clear a priori to which order in v the calculations must 
be taken in order to obtain results of sufficient accuracy. 
To answer this question requires a different method of 
analysis. Perturbation theory is ideally suited. 

In perturbation theory it is assumed that the binary's 
mass ratio is very small, ^/M <C 1, and that the massive 
companion is a black hole (rotating or nonrotating) . On 
the other hand, the orbital velocity is not restricted: the 
motion is arbitrarily fast, and the gravitational field arbi- 
trarily strong. To calculate, within perturbation theory, 
the gravitational waves emitted by such binary systems 
is the purpose of this series of papers ]2l| , |27| -p0t , as well 
as that of the work reviewed below. 

The perturbation-theory approach can be summarized 
as follows. 

A particle of mass /i moves in the gravitational field 
of a black hole of mass M . (By "particle" we mean 
an object whose internal structure is of no relevance 
to the problem.) The particle possesses a stress-energy 



The theory of black-hole perturbations (gravitational 
and otherwise) has been the topic of a vast literature. A 
self-contained summary can be found in Chandrasekhar's 
book [|l|. In this series of papers we have adopted 
Teukolsky's formulation of perturbation theory [Q, in 
which the perturbation field is ^"4, a complex- valued com- 
ponent of the Weyl tensor. The equation satisfied by 
is known as the Teukolsky equation. 

The perturbation formalism was summarized in paper 
I and applied to the specific case of circular motion 
around a nonrotating black hole. Methods were devel- 
oped for integrating the wave equation analytically in the 
low- frequency limit, and the energy loss was calculated to 
order v'^ in the post-Newtonian expansion. These analyt- 
ical methods were extended by Sasaki [Q, and Tagoshi 
and Sasaki ||l^ have obtained E to posf'-Newtonian or- 
der [0(w*); see Sec. V D]. Much was learned from this 
work about the structure of the post-Newtonian expan- 
sion. 

In paper II the perturbation equations were inte- 
grated numerically, and the energy loss was calculated ex- 
actly (apart from numerical error) for a wide range of or- 
bital frequencies. By fitting a post-Newtonian expansion 
to the curve E{f) the quality of the approximation could 
be determined. It was concluded that for the purpose of 
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matched filtering, the expansion must be extended so as 
to include terms of order or higher. These numerical 
calculations were repeated with much higher accuracy by 
Tagoshi and Nakamura who confirmed this conclu- 
sion. The analytical work of Tagoshi and Sasaki |]l8| is 
in complete agreement with the numerical work. 

The gravitational waves emitted by a particle in eccen- 
tric motion around a Schwarzschild black hole were calcu- 
lated, and the orbital evolution under radiation reaction 
considered, by Tanaka et al. [Q. They discovered that 
while radiation reaction always reduces the eccentricity 
in slow-motion situations, strong- field radiation reaction 
can cause the eccentricity to increase. This was studied 
in detail in paper III ||2^ for the special case of small 
eccentricities, and by Cutler, Kennefick and Poisson ||3C| ] 
(we shall treat this paper as a honorary member of this 
series) for the general case. This increase of eccentricity 
can be understood as a precursor effect to the eventual 
plunging of the orbit, at the end of the inspiral. 

Generalization of the work described above to the case 
of a rotating black hole was the topic of several papers. 
In the limit of slow rotation, analytical methods were 
used in paper IV |2|] to calculate the corrections to the 
wave forms and energy loss due to the black hole's rota- 
tion. These results were generalized by Shibata et al. , 
who considered a rapidly rotating black hole, as well as 
orbits slightly inclined with respect to the hole's equato- 
rial plane. Both papers assumed circular orbits. These 
results were extended to the case of slightly eccentric or- 
bits by Tagoshi |^ . Numerical results were obtained by 
Shibata |36|, who also considered the case of eccentric 
orbits 



E. Black-hole absorption 

We now turn to a description of the work contained in 
this paper. 

As mentioned previously, solving the perturbation 
equations amounts to solving a wave equation in the 
background of a Schwarzschild black hole. Boundary con- 
ditions must be imposed when integrating this equation. 
The correct choice is to impose a no-incoming-radiation 
condition, which forces the system to lose energy to, and 
not gain energy from, gravitational waves. In practice, 
this condition is implemented by imposing outgoing- wave 
boundary conditions at infinity, and ingoing- wave bound- 
ary conditions at the black-hole horizon. 

In the analytic work of Poisson , Sasaki , and 
Tagoshi and Sasaki |l8| , the horizon boundary conditions 
were not found to play a significant role. The black hole 
could have been replaced by some regular distribution 
of matter, and the horizon boundary conditions replaced 
by a regularity condition at the origin (this is the type 
of condition which is imposed in post-Newtonian calcu- 
lations (l^), and the results would have been identical. 
Due to this situation, the question of correctly imposing 



the horizon boundary conditions was not fully examined. 

The apparent irrelevance of the horizon boundary con- 
ditions turns out to be a consequence of two facts. The 
first is that only E°° , the energy radiated to infinity, was 
calculated in the previous analytical work. The total en- 
ergy loss, must also include ^ the energy flowing 
through the black-hole horizon. The second is that E°° 
was calculated only up to order in the post-Newtonian 
expansion. Sasak i |33l has shown [this is essentially a con- 
sequence of Eq. ( |1.4D below] that the horizon boundary 
conditions affect E^" only at order and beyond. 

One of the main objectives of this paper is to calculate 
E^ , the rate at which the black hole absorbs energy, 
to leading order in a slow-motion approximation. This 
calculation is carried out in Sec. V, using material derived 
in preceding sections. We find 
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so that contributes terms of order and higher to 
the post-Newtonian expansion of E. It can therefore be 
said that the black-hole absorption is an effect occurring 
at post^-Newtonian order in the energy loss. This effect 
utterly dominates that of the horizon boundary condi- 
tions on the energy radiated to infinity. 

Nevertheless, the black-hole absorption is a small ef- 
fect. In view of the requirement for matched filtering 
— energy loss accurate to post'^-Newtonian order ^9 27 
(see subsection D) — this effect can be safely ignored. 
The smallness of E"/E°° can be attributed to the pres- 
ence of a potential barrier in the vicinity of the black 
hole. This barrier is a manifestation of the spacetime 
curvature; it influences the propagation of waves near the 
black hole. Since gravitational waves generated by binary 
motion have a frequency / such that irMf = ^ 1, 
most of the waves propagating initially toward the black 
hole are reflected off the potential barrier (see subsection 
F). As a result, the black hole is unable to absorb much. 

That the black-hole absorption is a small effect is ulti- 
mately due to the extreme degree of compactness of black 
holes. How large the absorption would be for a (sub- 
stantially less compact) neutron star is unknown. This 
question warrants further examination. 



F. Integration of the Regge- Wheeler equation 

While integrating the perturbation equations amounts 
to solving a wave equation, solving that equation 
amounts to integrating the Regge- Wheeler equation 
p8[ for the radial function A"^£(r). Here, r is the 
Schwarzschild radial coordinate, w the wave's angular fre- 
quency, and i the spherical-harmonic index. The Regge- 
Wheeler equation takes the form 



V{r) 



(1.3) 
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where d/dr* = (1 - 2M/r)d/dr and V{r) = (1 - 
2M/r)[e{e + \)lr^ - GM/r^]. A significant portion of 
this paper (Sec. Ill) is d evoted to developing techniques 
for integrating Eq. (1^). These techniques rely on an 
assumption of low frequency, Mlo <C 1, which is appro- 
priate for slow-motion situations. 

Our techniques allow us to calculate the reflection and 
transmission coefficients associated with the potential 
barrier V{r). Suppose an incident wave e~''^(*+'' ' is 
sent from r = oo toward the black hole. The wave par- 
tially reflects off the potential barrier. The reflected wave 
^L^£e~'"*^*~'' ^ comes back to r = cx), while the trans- 
mitted wave 7^£e~'"(*"'"'' ^ goes through the black-hole 
horizon. In Sec. Ill E we find that the reflection and 
transmission coefficients are given by 



, - (phasc){l + 0[|Mw|2(^+i)] }, (1.4) 



and 



2(£- 2)1(^-1-2)! 
(2^)!(2^-hl)!! 

X |l + 'kM\u\ 4 



2Afw 1^+1 (phase) 
0[(A/c.)2]}. (1.5) 



Expressions for the phase factors, accurate to order A/w, 
can be found in Sec. Ill E. (In the language of that sec- 
tion, Tl^t = ^°"7A'", and = 1/A'".) 

Leading-order expressions for the reflection and trans- 
mission coefficients were previously derived by Fackerell 
|9|. (See also Price Q and Thornc |4lL) His results 
correspond to the Mw — > hmit of Eqs. (|1.4|) and (1.5), 
and therefore omit the OiMuS) part of the phase factors 
and the ttMIwI correction. These corrections, we should 



point out, are not needed for the derivation of Eq. (1.2) 



G. Perturbation theory as a wave generation 
formalism 

It is interesting to compare the wave generation formal- 
ism which derives from perturbation theory to the post- 
Newtonian formalism of Blanchet and Damour p"5|-p^ . 

The Teukolsky equation can be separated by de- 
composing the perturbation field ^'4 into spin-weighted 
spherical harmonics |^^. From the solution at large 
distances it is possible to reconstruct the (traceless- 
transverse) gravitational- wave field , which is then 
naturally expanded in tensor spherical harmonics (Sec. II 
B). Apart from an overall factor of l/r, the coefficients 
are functions of retarded time only, and play the role 
of radiative multipole moments. These moments can be 
expressed in terms of an integral over the source. In gen- 
eral the radiative moments are not simply related to the 
source moments. However, to leading order in u in a 
slow-motion approximation, the relation is simple. This 
is estabhshed in Sec. IV A. 

The situation is analogous in the Blanchet-Damour for- 
malism. Here, the radiative field is also expanded in 



multipole moments, but for convenience a symmetric- 
tracefree (STF) representation is favored. (An expansion 
in STF moments is entirely equivalent to an expansion in 
tensor spherical harmonics [[43[ .) In general the radiative 
moments are not simply related to the source moments, 
except for the limiting case of slowly moving sources. 

The Blanchet-Damour formalism uses two coordinate 
systems. The first, {t, a;}, is rooted to the source, and 
chosen so as to make near- zone calculations simple. How- 
ever, the true light cones of the spacetime do not every- 
where coincide with the near-zone light cones (described 
by i — |a;|), and this generates artificial divergences in the 
radiation field. To remedy this, Blanchet and Damour 
p^t introduce radiative coordinates, {T, X}, which are 
adapted to the true light-cone structure of the space- 
time. In terms of these coordinates, the radiation field is 
finite. The coordinate systems are approximately related 
by T- \X\ = t- \x\-2M\u(\x\ls), where effects nonlin- 
ear in the gravitational-wave field have been neglected. 
The log term is naturally interpreted as a Shapiro time 
delay, and s is an arbitrary parameter which serves to fix 
the origin of T . 

No such coordinate transformation is needed in per- 
turbation theory. One can work with the usual 
Schwarzschild coordinates, and the radiative multipole 
moments are naturally expressed as functions of retarded 
time u = t-r* = i - r - 2Mln(r/2M - 1); see Sec. II B. 
The Shapiro time delay is therefore automatically incor- 
porated into the retarded time, and the radiative field is 
automatically finite. 

The Blanchet-Damour formalism has been used to cal- 
culate the tail part of the radiation field |l^. Such 
tails arise because the curvature of spacetime outside the 
source scatters the gravitational waves as they propagate 
toward the far zone. As a result, the radiation field at 
time T does not depend only on the state of the source at 
retarded time T — \X\, but also on its entire past history. 
This effect modifies both the amplitude and the phase of 
the radiative multipole moments. 

The tail effect can also be investigated using pertur- 
bation theory, which we do in Sec. IV B. The results are 
identical to that of post-Newtonian theory. 

We therefore see that perturbation theory can be used 
to recognize and clarify many of the issues that must 
be addressed when dealing with post-Newtonian wave 
generation. This, we feel, is a useful aspect of our work. 



H. Organization of the paper 

We begin in Sec. II with a brief summary of the pertur- 
bation formalism. This formalism has been manipulated, 
in the course of this series of papers, so as to make it as 
easy to use as possible. We believe that the present for- 
mulation is as convenient as can possibly be. We will not 
present derivations in this section, but refer the reader 
to earlier references. Section II A is devoted to a de- 
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scription of the inhomogeneous Teukolsky equation and 
its integration by means of a Green's function. A proce- 
dure to obtain this Green's function in terms of solutions 
to the Regge- Wheeler equation is described in schematic 
terms. The procedure is detailed in Appendix A. Section 
II B describes how to obtain the gravitational-wave field 
from the large-distance behavior of the Teukolsky func- 
tion. The radiative multipole moments are introduced. 
Finally, Sec. II C gives expressions for the rates at which 
the gravitational waves carry energy and angular momen- 
tum, both to infinity and down the black hole. 

In Sec. Ill we present our techniques for integrating 
the Regge- Wheeler equation in the limit of low frequen- 
cies. Two sets of boundary conditions are imposed. The 
function X^^{r) satisfies ingoing- wave boundary condi- 
tions at the black- hole horizon, while X^^jr) satisfies 
outgoing- wave boundary conditions at infinity . Most 
of the section is devoted to X^^{r). Some notation is in- 
troduced in Sec. Ill A. The Regge- Wheeler equation is 
integrated near r — 2M in Sec. Ill B, and is integrated 
in the limit Muj ^ 1 in Sec. Ill C. Matching is carried 
out in Sec. Ill D. The reflection and transmission coeffi- 
cients are then calculated in Sec. Ill E. Finally, Sec. Ill 
F contains a brief discussion of the function Ar^(r). 

In Sec. IV we illustrate how black-hole perturbation 
theory can be used as a wave generation formalism. We 
consider slowly moving sources, and first (Sec. IV A) ob- 
tain expressions, valid to leading order in v, for the radia- 
tive multipole moments. Then, in Sec. IV B, we consider 
the corrections to the radiative moments which are due 
to wave-propagation (tail) effects. 

In Sec. V we present calculations pertaining to gravi- 
tational waves produced by a particle in circular motion 
around a black hole. These calculations are carried out 
within the slow-motion approximation, to leading order 
in V. In Sees. V A and B we derive expressions for the ra- 
diative multipole moments. In Sec. V C we calculate the 
contribution from each moment to the energy radiated. 
We consider both the flux at inflnity, and that at the 
black-hole horizon. Finally, in Sec. V D, we consider the 
role of the black-hole absorption in the post-Newtonian 
expansion of the energy loss. 



A. The Teukolsky equation and its solution 

In the Teukolsky formalism, gravitational pertur- 
bations of the Schwarzschild black hole are rep- 
resented by the complex-valued function \E'4 = 
—Cafi'y5'n°"m^n''fh^- Here, Ca/i-yS is the Weyl tensor, 
and = i(l,-/,0,0), = (0, 0, 1, -i csc6')/V2r 
(in the {t, r, 9, </>} Schwarzschild coordinates) are mem- 
bers of an orthonormal null tetrad. We have introduced 
/ = 1 — 2M/r, and a bar denotes complex conjugation. 

The Weyl scalar can be decomposed into Fourier- 
harmonic components according to 



1 f°° 



)e-''^*, (2.1) 



where sYemid-,4') s-^e spin- weighted spherical harmonics 
Q . The sums over £ and m are restricted to —£ < m < i 
and £ > 2. The radial function R^im{f) satisfies the 
inhomogeneous Teukolsky equation [R2| , 



^'/rf;:2-2(r-M)- + t/(r 



dr 



(2.2) 



with 

U{r) = f-^ [{ujrf - iiujir - 3M )] - {£ ~ l){£ + 2). 

(2.3) 



The source term to the right-hand side of Eq. (2.2) is 
constructed from the particle's stress-energy tensor, 



T"^(.t) =fi dr u'^u'^S'^^^ [x - x'{t)] , 



(2.4) 



where x represents the spacetime event and x'{t) the 
particle's world line with tangent vector u" = dx'^/dr 
(t denotes proper time). The first step is to obtain the 
projections qT — TafjU^n^ , -iT — Tapn'^rh^ , and -2T = 
Tapfh^fh^ . Then one calculates the Fourier- harmonic 
components sTuimir) according to 



II. THE PERTURBATION FORMALISM 



sTuiim{r) = — dt dO. sT sYi„i\ 

2,71 



(2.5) 



The stress-energy tensor associated with a moving par- 
ticle perturbs the gravitational field of a nonrotating 
black hole. We describe this perturbation using the 
Teukolsky formalism js^, which we review below. The 
following presentation will be brief; missing details can 
be found in Refs. pl|,p8|j30| , and in Appendix A. 



where dO. is the element of solid angle. Finally, T^imir) 
is obtained by applying a certain differential operator 
to each sTuJtm{f')^ and then summing over s pl| , ^ . 
Schematically, T — sDgT, where gD (other indices 
suppressed for notational simplicity) are the differential 

operators. (See Appendix A for details.) 

The inhomogeneous Teukolsky equation ( |2.2|) can be 
integrated by means of a Green's function [46|]. The 
Green's function is constructed from two linearly in- 
dependent solutions to the homogeneous equation, so 
that \E'4 satisfies a no-incoming-radiation condition. 
Schematically (indices suppressed), we have G{r,r') = 
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{r<^)R°°{ry), where G{r,r') is the Green's function, 
R^{r) the solution to the homogeneous equation sat- 
isfying ingoing-wave boundary conditions at the black- 
hole horizon, R°°{r) the solution satisfying outgoing- 
wave boundary conditions at infinity, and r< (r>) de- 
notes the lesser (greater) of r and r' . Schematically also, 
the solution takes the form R{r) = J dr' G{r, r')T{r') — 
j:jdr'G{r,r')sD,T{r'). 

It is then useful to define the adjoint operators 
so that R(r) can be more conveniently expressed as 
Rir) = Y.s S dr' sT{r')sD^G{r,r'). The last step con- 
sists of invoking the Chandrasekhar transformation , 
which relates solutions to the homogeneous Teukolsky 
equation to that of the Regge- Wheeler equation 
Thus, G{r, r') can be conveniently written in terms of lin- 
early independent solutions of the Regge- Wheeler equa- 
tion. (Relevant details can be found in Appendix A.) 

The Regge- Wheeler equation psl takes the form 



Vir) 



(2.6) 



where d/dr* = fd/dr and V{r) = f[i{t-l)/r'^-QM/r% 
The solutions of interest are Ar^(r), which satisfies 
ingoing-wave boundary conditions at the black-hole hori- 
zon, and X'^g (r) , which satisfies outgoing- wave boundary 
conditions at infinity |^] . More precisely, these functions 
are defined so as to have the following asymptotic behav- 
iors: 



X"e.{ 



X, 



(r 



2M) 

-> oo) 
oo) 



e 

A" 



A out Aujr* 



(2.7) 



Here, r* = r + 2A/ln(r/2Af - 1), and A]^^ and A°f are 
constants. It follows from the conservation of the Wron- 
skian that X^(r 2M) e"'"'^''' -I- Aj^"^e*"''*. 

Schematically, the Chandrasekhar transformation |^] 
takes the form R^'°° — CX^'°^, where C is a second- 
order differential operator. Thus, the effective Green's 
functions sD^G{r,r') can be expressed in terms of 
X^'°°{r) and their derivatives. Because the Regge- 
Wheeler functions satisfy a second-order differential 
equation, the final expressions involve the functions and 
their first derivatives only. (See Appendix A for details.) 

At large radii, the solution to the inhomogcncous 
Teukolsky equation is given by 



Ruimir oo) ^ ^^io^Z^,^r''e''^''' . 



Near the black-hole horizon |28 



Ru. 



(2.8) 



(2.9) 



After carrying out the manipulations described above. 



the "amplitudes" Z^^^ are found to be given by 



ryH.OO 



H.oo /I in 



X / drrf-\T^e^ir),r^eX^i°^{r). (2.10) 

We have introduced the following symbols: 

'^"i - 3 [(^ - +!)(£ + 2) - 12zMc.] , 



(2.11) 



-16(1 - 2iMuj){l - AiMuj){Mujf 



and 



2[{i-lW+l){li + 2)\'''^ s = 

sPi= \ 2[2{t-l){i + 2)\^''^ s = -l 

1 s = -2 

We also have 



(2.12) 



0^ Ull 



2(1 - Ml/r + iujr)rf— + f UU + 1) - QM/r\ 
dr ■' 

+ 2iLor{l - 3M/r + iujr), 
-~f {[£{£+ l) + 2iujr]rf 



dr 



^ulf. — 



-2-1- o) 



+ £{£ + l){f + iujr) - 2iujrf 
fl^2[{£-l){£ + 2) + 6M/r]rf 
+ {£- 1){£ + 2) [£{£ + 1) + 2iujr] + 12/M/r 



(2.13) 



d_ 

dr 



As noted previously, s^uje are first-order differential op- 
erators. That this is so is convenient both for analytical 
and numerical evaluation of Z, 



Ha 



B. Waveform and multipole moments 

The gravitational-wave field can be obtained from the 
asymptotic behavior of ^"4 at large distances. Choosing 
the 9 and (p directions as polarization axes, the two fun- 
damental polarizations of the gravitational waves can be 
expressed as [ESl 



h = — ih; 



2/Lt 



where u ^ t — r* , and 



dujZ"e-'^". 



(2.14) 



(2.15) 



The traceless-transverse gravitational- wave tensor is then 
hj^ = hmamb + hrharfib, (2-16) 
where latin indices denote spatial components. 
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It is clear that the Z^{u) represent the radiative mul- 
tipole moments of the gravitational field. To relate these 
quantities to t he com monl y used definitions, we substi- 
tute Eq. ( 2.14 ) into ( 2.f 6 ) and rewrite in terms of the 
Mathews tensor spherical harmonics |^^. The result is 
Eq. (4.3) of Ref. ||, 



. TT 



- I \rT^E2.ln 



■ Sp. 



( \rpB11-n 



(2.17) 



where the T-tensors are the spherical harmonics. We 
thus find that the mass multiple moments are given by 

X,^(u) = yf2pi\^f^{u) + (-f)'"I^,,_„(^i)], (2.18) 

while the current multipole moments are given by 

5ft„(M) = ^^2^l\^Zl^(u) - (-f)"I^,,_„(?/)] . (2.19) 

When the source is confined to the equatorial plane, 
the identity ^ = (~l)^^(5fm holds (see paper III 
psf for a derivation). This simplifies the expressions for 
the mass and current multipole moments. 



III. INTEGRATION OF THE REGGE- WHEELER 
EQUATION 



Our objective in this section is to integrate Eq. ( |2.6| ) 
in the limit Mlo <C 1. Techniques for solving this prob- 
lem were developed by Poisson in the first paper in this 
series ||2^] , and then extended by Sasaki . We extend 
these techniques further here, by explicitly considering 
the issue of the boundary conditions at the black-hole 
horizon. 

Most of this section will be devoted to the "ingoing- 
wave" Regge- Wheeler function X^(r). The "outgoing- 
wave" function X'^gir) will be briefly considered in sub- 
section F. For notational simplicity we will, throughout 
this section, suppress the use of the ujI suffix. 



A. Preliminaries 

For concreteness we assume that the wave frequency 
Lo is positive. Results for negative frequencies must be 
derived separately |pO| . For convenience we define 



C. Energy and angular momentum fluxes 



Equations (2_^) and ( p.9| ) can be used to calculate the 
rates at which energy and angular momentum are radi- 
ated to infinity and absorbed by the black hole. Here 
we consider only the special case of circular orbits; more 
general results can be found in Ref. p^ . 

It will be shown in Sec. V that for circular orbits, the 
frequency spectrum of the gravitational waves contains 
all the harmonics of the orbital frequency = d4)/dt, so 
that 



jH.oo 



7H,, 



^S{u! — mQ). 



(2.20) 



For this special case, the energy flux at infinity is given 
by @ 



E°° ^ 



i(^)'E(^M^I^fj^ (2.21) 



while the energy flux at the event horizon is |48| 



An \M 



Im 



at-n 



-'iml 



where 



2i2[l + 4(Mcj)2] [l + 16(Afa;)2] 
I {i - l)i{e + l)ie + 2) - 12iMLu\ 



(2.22) 



[Mujf. (2.23) 



In Eqs. ( 2.21 ) and ( 2.22 ), a dot denotes differentiation 
with respect to coordinate time t, and oj = mfl. The total 
amount of energy radiated by the gravitational waves is 
E = E'^ +E". 

The angular momentum fluxes can be obtained from 
the relations E°° = nL°° and E" = nL" , vahd for any 
radiating system in rigid rotation Q . 



z = LOT, £ — 2Mlu, z* = z + e ln(z — e) 



(3.1) 



Both z and e are dimcnsionlcss, and we n ote that z* — 
LOT* -\- elne; r* was introduced in Eq. (2.7). 

With these definitions the Regge- Wheeler equation be- 
comes 



dz*^ 



1-/ 



£(^-1-1) 3e 



Xiz) = 0, (3.2) 



where / = 1 — e/z. The ingoing- wave Regge- Wheeler 
function satisfies the following boundary conditions: 



X"{z^e) 



(3.3) 



To carry out the integration of the Regge- Wheeler 
equation it is useful to follow Sasaki [|3| and introduce 
an auxiliary function (z) defined by 



Y"{z) = z"ie*^'"("-'^)X^(z). 



(3.4) 



The great advantage of deahng with Y^{z) instead of 
(z) directly will be become apparent below. The aux- 
iliary function satisfies 



z{z 



£)— + [2{l - ie)z - e\ 



dz"^ 



dz 



z^ - £(£+!) + -(A- iz + z') '>Y"(z)^0. (3.5) 

z 



Equations (3.3) and 



further imply 



(3.6) 
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and 



j[in^ie Ine 



A 



outg — In e ^ — l^i{z-\-2e In z) 



(3.7) 



For future reference we shall now rewrite Eq. (3/?) in 
a different form. We first expand the right-hand side in 
powers of e, which below will treated as a small number. 
(This assumption has not yet been used.) We also invoke 
the spherical Bessel functions jt{z) and (z), which will 
play a prominent role below, via the asymptotic relations 



obtain 

Y"{z- 



(±z) + [jf(z — > oo) ± int{z — > oo)]. We thus 



[a+ + 2melnz + 0{e'^)\ji{z —> oo) 
— e[ia_ + 2alnz + 0{e)\ni{z oo). 



In terms of a± we have 
A"' = i(i)'+i(a 



-ie \n€ 



ea_je 

= i(-i)^+i(a+ - ea_)e'^'"^ 
a = h{a+ — ea_). 



(3.8) 



(3.9) 



B. Auxiliary function: solution for z <^ 1 



We first integrate Eq. ( B.f: ) in the limit z <C 1 . We shall 
also treat e as a small number, but leave arbitrary the 
ratio z/e. Accordingly, we neg lect the terms within 
the large square brackets of Eq. ( |3.5| ): the first one can be 
neglected in front oi £{£+!), while the second is negligible 
in front of 4. All other terms are kept in the equation. 

The resulting equation can easily be solved if we change 
the independent variable to a; = 1 — z/e and the depen- 
dent variable to Z = (e/z)^F. We then obtain 

x{x - l)Z" + [2(3 - ie)x - (I - 2ie)] Z' 

+ [6-i{£+l)-5ie]Z = 0, (3.10) 

where a prime denotes differentiation with respect to x. 
This is the hypergeometric equation, with parameters 



b = i + 3-ie + 0{e^), 
c = 1 - 2is. 



(3.11) 



The two linearly independent solutions are F{a,b]c;x) 
and x^~'^F{a + 1 — c,b + 1 — c;2 — c; x), where F is the 
hypergeometric function. However, only the first solu- 
tion is regular at x = {z = e). Compatibility with 
Eq. (^^) therefore implies that the second solution must 
be rejected. We therefore obtain 

y^(z < 1) = e-ie*^('""-iHz/£)2F(a, 6; c; 1 - z/e). 

(3.12) 



Equation (3.12) enforces the correct boundary condition 
at the black-hole horizon. 

For later use we now evaluate Y^{z) in the limit 
£ ^ z ^ 1. First, a straightforward calculation, us- 
ing Eg. (1 5.3.7) of Abramowitz and Stegun reduces 
Eq. (|3.12|) to 



Y"{e < z < 1) e 



X 1 



r(c)r(6-a) /zxWs 
r(6)r(c-a) Ve. 
(^-2)(^ + 2) 



2£ 



le 



0[{e/z) 



(3.13) 



Next, we expand the complex exponential, and the F 
functions according to F(?i ~\- 5) — (n — 1)![1 -I- V'("-)'5 + 
0((5^)], in powers of e. [Here, n is an integer and ■(/'(n) = 
—"f+Y^=i 1 where 7 ~ 0.57721 is the Euler constant, 
is the digamma function.] Finally, we obtain 



y^(£<z< 1) 



{2£)\ 



(^-2)!(^ + 2)!e^+i 
{l-2){l + 2)_e_^ 
z 



1 + ie{ai + In z) 



21 



(3.14) 



We have defined 

at = 2-1 + - 1) -t- + 3) - 1. 

C. Auxiliary function: solution for e ^ 1 



(3.15) 



We now integrate Eq. (3.5) in the limit e ^ 1. It will 



also be assumed that z ^ e, but z is no longer restricted 
to be much smaller than unity. 

We first rewrite Eq. (3.5) by sending to the right-hand 
side all terms which are linear in e; the left-hand side is 
then independent of e. To solve this equation we proceed 
iteratively Q, by setting 



(3.16) 



A short calculation then shows that each (z) satis- 

fies 



dz^ 



+ 2z— + z"^ - £{t+l) 



dz 



r^(")(^) = ^^'"Hz), 

(3.17) 



where 

5^(")(z) 



z-^^ + z(l + 2iz) 



dz2 

{A~iz + z^) 



d 
dz 



y^("-i)(z). (3.18) 
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Equation (3.17) is an inhomogeneous spherical Bessel 
equation. It is the simphcity of this equation which mo- 
tivated the introduction of the auxihary function |33| . 

The zeroth-order solution, Y^^^'>{z), satisfies the ho- 
mogeneous spherical Bessel equation, and must therefore 
be a linear com bination of je{z) and ni[z). Compatibility 
with Eq. ( ^.14 ) then demands 



(3.19) 



where B is constant to be determined. 

The procedure to obtain F^'^Hz) was described in de- 
tail in Ref. |3^. In short, Eq. (3.19) is substituted into 
Eq. (|3.18|) to obtain S"'^^\ and Eq. (|3.17|) is integrated 
using the Green's function G{z,z') = j£(z<)n£(z>). The 
resulting integrals can be evaluated explicitly, and the 
final result is 

Y"^^\z)/B^ (a^i) +ilnz-Si2z)jf 

{h^^^ -hi2z + G\2z^ni 

-2 



+ - + — — r z^{ndp - .nnp)jp 



21-1 



{l-2){l + 2) 

2l{2t+l) '^£{£-1) 



Je-i 



2(£-M)(2£-|-l) 



(3.20) 



Here, a^^^ and b^^^ are constants of integration to be de- 
termined. Six = dtt~^ sint is the sine integral, and 
Cix = 7-|-lna;-|- dtt^^{cost — 1) is the cosine integral. 

D. Matching 



The constants B, a^^' and 6'^^-' can be determined by 
matching y^(z) = y^(o)(z) + eY"^^\z) + O(e^) to 
Eq. (13.14 ). This involves the evaluation of Eqs. ( |3.19| ) 
and (3.20) in the limit z <C 1. This calculation is straight- 
forward, since the limiting expressions for the spherical 
Bessel functions, as well as that for the sine and cosine 
integrals, are well known. We find that — —7 and 
that a*^^' is arbitrary — it can be absorbed into the 0{e) 
part of B. Without loss of generality we may set a*^^^ = 0. 

We no te in passi ng that consistency between 
Eqs. ( |3.14| ), ( |3.19| ) and ( |3.20| ) demands that the following 
equation hold: 



CXJ ^ 

■^-^ 2n(2n 

n=2 ^ ' 



1-2 . 



1 



Z2'[l + 0(z2 



2 • 2 

Z Jp 



(2^- l)!!(2^-3)!!^(£- 1)' 



(3.21) 



We hav e not been able to directly establish the validity 
of Eq. ( ^.21 ), but have checked that it does indeed hold 
for several special cases. 

Finally, matching yields the value of B, which is 



B 



(2^)!(2^-h 1)!! 1 
{£ - 2)\{e + 2)l7^ 



j[l + ieae + 0{e^)], (3.22) 



where ai was introduced in Eq. (3.15). The integration of 
Eq. ( |3.5D , accurately to first order in e and with boundary 
condition ( ), is now completed. In Ref. |Q , Sasaki has 
indicated how to proceed to higher order in s. 



E. Calculation of A"" and 

We are now in a position to evalua te Y ^(z) in the 
limit z —^ 00 and to compare with Eq. (3.8) so as to ob- 
tain expressions for the amplitudes A™ and Again 
the calculation is straightforward, since the asymptotic 
expressions for the spherical Bessel functions (as well as 
that for the sine and cosine integrals) are well known. 

We use such expressions as ji^i ~ —je+i ~ —ne and 
z'^inijp - jinp)jp ~ i[l - {-lY~P]ne to calculate 



Y"{z 



B 



{[l + e{tlnz-TT/2)]ji{z 
e(ln2z — Pi)ne{z 00) 



(3.23) 



where 



i^{£) + ^{e + i) 



li£ + l) 



(3.24) 



Comparison with Eqs. (3.8) and (|3.9|) finally yields 



A" 



(2£)!(2£+l)!! fi 
2{£ - 2y.{£ + 2)1 (1 



i+1 



-is (In 2e — ai — 



(3.25) 



^out . 



(2^)1(2^+1) 
2(^-2)!(^-H2)! 



M /-i 



i+1 



1 ) iE{ln2e+ae-pe) 



Analytical expressions for A™ and in the low- 

frequency regime were first obtained by Fack erell 39 1 ; his 
results correspond to the e ^ limit of Eq. ( 3.25 ). That 



these quantities scale like £~(^+i) in that regime was first 
discovered by Price and Thorne |^l|. The 0(e) cor- 
rections to Fackerell's results, we believe, are presented 
here for the first time. 

It is important to notice that ^4°"* is not quite the 
complex conjugate of A™. (The physical relevance of the 
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phase facto r wil l be discussed in Sec. IV B.) In magni- 
tude, Eq. ( |3.25| ) implies that they are equal, |^°"*| = 
+ 0(e^)], up to a fractional accuracy of order e^. 
In fact, a much stronger result follows from the Wron- 
skian relation — |v4°^*p = 1. Simple manipulations 

and use of Eq. (|3.25|) indeed reveal that |^°"*| and 
are equal up to a fractional accuracy of order e'^^^'^^\ 



eventually obtain that c^^'^ is arbitrary and can be set to 
zero without loss of generality, that = 7r/2, and that 



C = iiY+^ [1 + ze(/3, + 7 - Ine) + Oie^)] , 



(3.30) 



where f3e is given in Eq. ( 3.24 ). This completes the de- 
termination of y°°(z) to first order in e. 



F. The outgoing- wave Regge- Wheeler function 

We now integrate the Regge- Wheeler equation for 
X^^{r), the out going- wave function. The method is en- 
tirely analogous to that described in the preceding sub- 
sections. We briefly sketch the most important steps. 

We first introduce an auxiliary function Y°°{z) defined 

by 



y°°(z) 



^X°°iz). 



(3.26) 



From Eq. (2.7) we see that the auxiliary function must 
satisfy the boundary condition 



y°°(z^oo) 



(3.27) 



It is immediate that Y°°{z) satisfies the same differential 
equation, Eq. (3^), as Y^{z). Since the boundary con- 
ditio ns a re imposed at z = oo, we shall not need to solve 
Eq. ( |3.5| ) in the Hmit z ^ 1 (as described in subsection 
B). Instead we can proceed with an analysis similar to 
that contained in subsection C. 

The zeroth-order solution must be identified, up 

to a normalization constant, with the spherical Hankel 
function h^^"^ (z), whose asymptotic behavior is identical 



to Eq. (3.27). We therefore have 



y°°(")(z) = C/if'(z), 



(3.28) 



where C is a constant to be determined. Taking the com- 
plex conjugate, we get Y°°^^\z) = Ch^p{z) = C[je{z) — 
ini{z)]. The integration of Eq. (|3.17|) for y°°(i)(z) pro- 



ceeds as previously described. We obtain 

(z)/C = c(i'4^^ + - Si2z + iCi2z ) h 



,(2) 



P p+l 



p=l 



2e- 1 



2^(2^ + 1) iie-i) 

ji -!){£ + 3) (1) 
2(£+l)(2£+l) ^+1' 



(3.29) 



where c^^^ and d*^^^ ar e con stants of integration. 

Matching to Eq. (3.27) involves the evaluation of 
Eq. (3.29) in the limit z — s- oo, which requi res m anip- 
ulations similar to the ones leading to Eq. ( |I23| ). We 



IV. RADIATIVE MULTIPOLE MOMENTS IN 
THE SLOW-MOTION APPROXIMATION 

In this section we apply the results of Sees. II and III 
to the limiting case of a slowly moving source. We first 
derive leading-order (Newtonian) expressions for the ra- 
diative multipole moments — Z^^^ in the frequency do- 
main, Z^{u) in the time domain; see Sec. II B. We then 
improve on these expressions by considering the correc- 
tions due to wave-propagation (tail) effects, which are 
of post'^/2_]\je^tonian order. These corrections do not 
depend on the specific details of the source, but are a 
manifestation of the curvature of spacetime outside the 
source. The larger, source-specific, post-Newtonian cor- 
rections will not be considered in this paper. 

For the purpose of the following calculation we do not 
assume a specific form for the stress-energy tensor T"^ 
that perturbs the black hole's gravitational field. In par- 
ticular, we do not assume that the perturbations are pro- 
duced by a moving particle and that T"^ is of the form 
( ^.4|) . [The following results will nevertheless contain fac- 
tors of the inverse of the particle's mass. This is 
because such a factor has been inserted for convenience 



into the definition of Z, ; see Eq. ( pJO| ). Such meaning- 
less occurrences of can be avoided by setting /i = 1 
throughout this section.] We do assume, however, that 
the source is slowly moving, in a sense made precise be- 
low. 



A. Leading-order calculation 

Let the source be characterized by some internal ve- 
locity V. We demand f <C 1. We assume that the 
source is gravitationally bounded to the black hole, so 
that M/r « w^, where r is a typical value of the radial 
coordinate inside the source. Finally, if w is a typical 
frequency of the gravitational waves, then slow motion 
implies « w, and Muj « u'^. In terms of the variables 
introduced in Sec. Ill, we have z R 



', s/z ~_v\ and 
This implies that such results as Eqs. ( 3.14 ) and 



( |3.25 ) can be used. 

In the slow-motion approximation, the stress-energy 
tensor is dominated by the component qT — Tapn"'n^ = 
\p+0{pv), where p is the source's mass density. By com- 
parison, _iT = 0{pv) and -2T = 0(pv^)-, these com- 
ponents will be neglected. Equation (|2.5| ) then implies 
oTujimir) = jj dVtp[u}, x)Y£rniO, (f) + 0{pv), where 
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(4.1) 



is the mass density in the frequency domain. 

Th e slo w-moti on a pproximation also implies that 
Eqs. ( I2I1I ) and ( |2l^ ) reduce to ^ff^ = - + 
1)(£ + 2) + 0(v^) and oO = 2rd/dr + £(£ + 1) + 0{v). 
Equations (|3J), ( |3.14| ), and ( |3.25| ) can then be combined 
to give 

_ 2(f + i)(£ + 2) 

Aj^-, " (2^+1)!! ^ ^ ' ^ ' 

to leading order in v. Substituting these results into 
Eq. (I2.IOD yields 



7" - 



47r 



(2£+l)!!/i 



{£ + !){£ + 2) 



-,1/2 



{£-!)£ 



(4.3) 



the radiative moments in the frequency domain. 

Finally, taking the Fourier transform of Eq. (4.3) we 
obtain, to leading order in v, 



47r 



(2£+l)!!^ 



(l^\){l + 2)Y'^ ( d_ 
du 



X / d'^x p{u,xy Yi,niO,< 



(4.4) 



the radiative moments in the (retarded) time domain. As 
expected, the radiative moments of order £ are given, up 
to a numerical factor, by the ^-th time derivative of the 
source moments, / d'^x pr^Yim- This leading-order calcu- 
lation agrees with that of Sec. V C of Ref. p3] . (It should 
be noted that we only have obtained the mass multipole 
moments. To calculate the current moments would in- 
volve keeping terms of order pv in our approximations.) 



B. Tail correction 



given m Eq. ( p^ . This step is dictated by Eq. ( p4| ); 
according to our rules we may ignore the other correction 
terms, (iii) We modify by multiplying the leading- 
order expression by (1 - 7rM|w|)e-2^^'^(i"4^l"J-"*-'5') 



as dictated by Eq. ( |3.25D ; (ii is defined in Eq. ( |3.24| ). 

After simple manipulations, the final result is that the 
tail corrections to the radiative multipole moments take 
the form 

+ TrAf I^^De^^^^'^^'"''^^!'^!-^'). (4.5) 



Here, Z^^^ stands for the leading-order expression for 
the radiative moments. We have introduced pi = Pi — 
6[{£-l)£{£+l){£ + 2)]-\ or 



pe^^{£-l) + 



£^ + 7£^ + 12^-1-8 



(4.6) 



2^(^+ l)(^-h2) ' 

after using Eq. ( ^.24| ). 

A result identical to Eq. (^^) was first derived for 
the frequency-domain quadrupole {£ = 2) moments by 
Blanchet and Shafer |^^, and then gen eralized to mo- 
ments of arbitrary order by Blanchet 152]. Those results 
were derived within the context of post-Newtonian the- 
ory (see Sees. I B and G). Our results, which confirm that 
of Blanchet and Shafer, were obtained using a completely 
different method of analysis, that of black- hole perturba- 
tion theory. It is, of course, pleasing that such different 
methods yield identical results. 

Blanchet 's tail correction is identical to Eq. ( |4.5| ) 
except for a different representation of the multipole 
moments (he uses a symmetric-tracefree representation) , 
and a different phase. Blanchet's expression substitutes 
a constant K£ — 7 in place of p^, where p7|| 



e-2 ^ 
51 1 



fe=i 



2£2 -1-5^ + 4 
£{£ + !){£ + 2)' 



(4.7) 



It is easy to see that the two constants are related by 
pi = K£ — 7 -f 1/2. Our result therefore differs from 
Blanchet's by a phase factor e~**^'^ independent of £. 
This overall phase has no physical significance: it can be 
absorbed into a shift in the origin of the retarded time u. 



We now improve upon the leadin g-or der calculation 
by considering the corrections to Eq. (^^) which are due 
to wave-propagation effects (tails; for the purpose of this 
discussion it is simpler to work in the frequency dom ain) . 
More precisely, we consider the corrections to Eq. (4.5) 
which are of order Muj « v^, the smallness parameter 
appearing in the Regge- Wheeler equation. We shall ig- 
nore, for reasons given above, the corrections of order v 
and v"^. 

The steps are simple. We repeat the calculation de- 
scribed in the previous subsection, but with the following 
changes, (i) We rewrite Eq. ( p.llD as Acff^ = j{£ -!)£{£ + 
l)(£-h2)e-i2^*^'^/('^-i)^(^+i)(^+2), to a fractional accuracy 
of order (Mw)^. (ii) We similarly modify X^Ar) by mul- 
tiplying the leading-order expression by e^*"*-'*^'^, with ai 



V. CIRCULAR ORBITS IN THE SLOW-MOTION 
APPROXIMATION 

In this section we use the perturbation formalism to 
calculate the energy radiated by a particle in circular mo- 
tion around a Schwarzschild black hole. We perform this 
calculation in the slow-motion approximation, to leading 
order in the velocity v. We consider both the energy ra- 
diated to infinity, and that absorbed by the black hole. 
Part of the calculations presented in this section were also 
carried out in paper I (21|. We repeat these calculations 
here for completeness. The remaining calculations are 
presented here for the first time. Equations (5.17) and 
( p3.18 ) below were quoted without derivation in a footnote 
of paper III p8|. 
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A. Z^f^ for £ + m even 



We begin with the calculation of the amplitudes Z^i^a 
and to leading order in v (to be defined precisely 

below), for the case where ^ + m is an even integer. The 
case £ + m odd will be treated in subsection B. That 
these cases must be considered separately will become 
apparent below. 

As was discussed in Sec. IV A, the dominant compo- 
nent of the stress-energy tensor in a sl ow- motion approx- 
imation is i^T = Tapn"n^ . Equation (2.4) then yields 



qT = — d{r - ro)d(cos 



nt). (5.1) 



Here, tq is the radius of the circular orbit which, without 
loss of generality, has been put in the equatorial plane 9 = 
7r/2; u" is the four-velocity, and = d(j)/dt = (M/ro^)^/^ 
is the angular velocity. 
We define v as 



(5.2) 



and demand 
M/vq = V 
ro)6{cose)S{(j)-nt) 
tegration yields 



^ 1. Up to corrections of order 



^ Eq. (O) reduces to qT = (/^/4ro^)(5(r 



Substitution into Eq. (|2.5D and in- 



JJ_ 

4ro^ 



oY,m{^,0)S{r-ro)6{uj-mn). (5.3) 



That sTujem, and hence Z^p"^, is propo rtional to Siu — 
mVl) has been anticipated in Eq. ( ^.20| ). (This result is 
exact, and not a consequence of the slow-motion approx- 
imation. See Ref. ||2l[ for details.) We shall now factor 
out this 5 func tion, and work with the quantities 
defined in Eq. (|2.20D. 



Substituting Eqs. (fl.2D and (^) into ( |2.10D and inte- 
grating, we obtain 



Zf^ = 47r(-*) 



(2£+ 1)!! 
X oi"^m(f,0)«'. 



{£+l){i + 2) 



1/2 



(5.4) 



Equation ( p.4|) is valid up to fractional corrections of or- 
der (see Refs. for details.) 

We first use the re- 



A similar calculation yields Z^^. 



suits of Sec. Ill E to obtain X^(r) 

where Irip is the spherical Hankel function, to leading 
order in Muj = 0(v'\ Then, up to fractional cor- 
rections of order (wr)^ = 0(i;^), we have ^'^((j') = 
{iY{2£ — l)ll{u!r)~^ . Next, we use the le ading- order ex- 



pressions for oT^^e and A"'^, Eqs. (2.13) and (3.25), to 



obtain 



o^^jX^Hr) ^ o^+2, (e - mi - 2y.{£ + 2)! {MluY+^ 



A 



{2l)\{2l+l) 



[ujrY 



(5.5) 



to leading orde r in v. Finally, substitution of Eqs. (5.3) 
and dU) into (|2.10|) yields 



nfi(2l)\(2£^\) 
x{£- 2)l(£ + 2)\oYe^{^,0)v^'-\ (5.6) 

to leading order in v. 

An explicit expression for o^fm(f,0) can be obtained 
by expressing the spherical harmonics in terms of the 
associated Legendre polynomials, and using the fact that 
P™(0) = (-l)5(''-™)(£ + m- l)!!/(£ - m)!! if ^ -I- m is 
even, and Pp {0) = if ^ -f- m is odd pi . Thus, 



1/2 



Ol£™(f,0) = (-1)^(^+'") 



2^+1 
47r 



1/2 



,1/2 



[{£ - m)\{£ + my. 
{£ - my.\{£ + my.\ 



(5.7) 



is even, and oYem{^, 0) = otherwise. We there- 



if £ 

fore see that the dominant, s = 0, contribution to Z, 
vanishes identically if £ -|- m is odd. 



B. Z";^ tor £ + m odd 



H.oo 



If £ 



is odd the dominant contribution to Z, 



H,oo 
Im 



comes from the component _iT = Tcif:in°''fh^ of the 
stress-energy tensor. A calculation analogous to that of 
the previous subsection yields -iT = {ifiv/2\/2ro'^)S{r — 
ro)S{cos9)6{<j) — ilt), so that 



IflV 

2V2^ 



iYeU^,0)S{r~ro)S{uj-mn). 

(5.8) 



From Eq. ( |2.13| ) we have -iF^^ = -£{£ + l){rd/dr + 1). 
It follows that 



-^i-^f^'^'m^^i-rY^' (5.9) 



and 



_iF^,X-(r) 



^P^^, {£ -!)£{£+ I) 
' {2£y.{2£+l) 
X (Mw)^+i(cjr)-^. 



(£-2)1(^-1-2)! 

(5.10) 



Here, X^(r), X^(r), and were approximated as in 
the previous subsection. 
We arrive at 



Z^ 



H 



87r(-i) 



£ + 2 



£~l 



1/2 



{2£+iy.\ 



X _iy,„,(§,o)w 



(5.11) 
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and 



1^-2 ■ 



1/2 



m3(2£)!(2£- 



1) 



2)! 

(5.12) 



to leading order in v. 

An explicit expression for ^lYimi'^, 0) can be obtained 
as follows. First, we calculate the s = — 1 spherical har- 
monics II using = l)]-i/2(5e - 
icsc9d^)oYi„i{9,(j)). We then express the s = har- 
monics in terms of the associated Legendre polynomi- 
als. Since P™(0) = when -|- m is odd, we obtain 
-ii^f™(f ,0) - (-l)™[(2£+l)/47r€(£+l)]i/2[(f-m)!/(^-f 

m)!]^/^P™'(0) if i + m is odd. (Here, a prime denotes dif- 
ferentiation with respect to the argument.) Finally, use 
of the relation g P]^'{Q) = \P^+^{Q) - \{t + m)(£ - 
m-f l)P;"-\0) yields 



-il£™(f,0) = (-l)*(^+™-i) 



2^+1 

47r^(f + 1) 

{l-m)\\{t + m)\\ 
172 



1/2 



(5.13) 



if £ -f m is odd. For i + m even _iy£m(-|, 0) is nonzero, 
but its expression will not be needed. 



C. Energy radiated 

It is now a straightforward task to calculate the rate 
at which the gravitational waves remove energy from the 
system. We use the results of the preceding subsections, 
together with Eqs. ( 2.21| )-( p.23 ). For convenience we 
shall normalize our expressions for E"^ and E" to the 
quadrupole-formula expression, Eqy = (32/5)(/i/M)^u^'^ 



We therefore define the numbers 



H,oo 



such that 



j^H ,00 



32 / fi 



.10 \ " lH,QO 



Z-^ 2 



(5.14) 



The factor of 1/2 is inserted for convenience because of 
the symmetry ry^ 1^ = 77^,^°° which can be derived from 
= i-'^yz^e^ (see paper III g for a deriva- 
tion). 

Calculation yields, to leading order in v, 



T(2^+l)!!2 [oy^rni^^O)] V ^ 



(5.15) 



if -|- m is even, and 



57r 



[1 + 2) 



(i^Pif[-.^.™(f,0)]V(-^) (5.16) 



if ^ -|- m is odd. 

Calculation also yields, to leading order in w. 



Tit 



H 



2l\ 



2(1+2) 00 
^ Vim 



if £ + TO is even, and 



if - 



is odd. 



2(^+1)! 



2(1+2) oc 



(5.17) 



(5.18) 



D. Black- hole absorption 

In a slow-motion approximation, the dominant con- 
tribution to the energy radiated comes from the i — 2, 
\m\ = 2 (mass quadrupole) terms in Eq. (5.14). The dom- 
inant contribution to the black-hole absorption is given 
by '722 = ^^'722 = ■ Therefore 



E" /E°^ =i;«[l + 0(i;2)]. 



(5.19) 



(That the correction must be of order can be estab- 
lished by direct calculation.) 

In a slow-motion approximation, the total energy ra- 
diated, E — E°° + E^ , can be expressed as a post- 
Newtonian expansion of the form 



E = Ei 



QF 1 1 + Oiv^) + 0{v^) + Oiv^) + Oiv^) 
+ 0{v^) + 0{v^ \nv) + 0{v'^) + 0{v^) 
+ 0{v^\iiv) + ■■■]. (5.20) 



The first three terms of this expansion, through 0{v^), 
were calculated in paper I all other terms were cal- 
culated by Tagoshi and Sasaki [Q. (See Sec. I B and D 
for a more detailed discussion.) 

We therefore see that the contribution E^ to the to- 
tal energy radiated occurs at quite a high order in the 
post-Newtonian expansion: 0(v^), or posf^-Newtonian 
order. What is more, the black-hole absorption is also a 
small contribution to the term of order v^. The coeffi- 
cient of the 0{v^) term in the post-Newtonian expansion 
of E°° was calculated, in the limit of small mass ratios, 
by Tagoshi and Sasaki . They find that it is approxi- 
mately equal to -117.5044. The coefficient o f the 0{v^) 
term in the expansion of E^ is given in Eq. ( 5.19| ): it is 
equal to unity. The black-hole absorption therefore con- 



tributes less than one percent of the term of order v in 



Eq. 

The black-hole absorption is a small effect indeed. 
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APPENDIX A: DERIVATION OF EQ. (2.10) 

The source term to the right-hand side of the inhomo- 
geneous Teukolsky equation, Eq. (2^), is given explicitly 

by [ill 



(Al) 



where the constants sPe are listed in Eq. ( 2.12| ), the func- 
tions sTuiemif) given in Eq. (2.5), and where 



.s 



3 f-1 



(A2) 



rfCr f Cr s 



fd/dr + iu! 



are differential operators. Here, £ 
d/dr* -\- iLu. 

Equation ( ^.^ ) is integrated by constructing a Green's 
function from two linearly independent solutions to the 
homogeneous equation |46|. These are denoted R^({r) 
and R'^({r), and have the following asymptotic behav- 
iors: R^,{r -> 2M) ^ (^r)V'e— i?5,(r -> oo) ^ 
Q'^^{ujr)-^e-'^'-' + QZfiujrfe'^"-' (Q'^^ and Q^f are 
constants); R'^f{r oo) ~ (cjr)^e'"''' . A straightfor- 
ward application of the general theory of Green's func- 
tions then shows that the solution to the inho mog eneous 
Teukolsky equation reduces to Eqs. (|2.8| ) and (2^), with 



drr-V-'CrWr.^mM- (A3) 



2M 



Our go al in this Appendix is to rewrite Eq. (|A3D into 
the form ( |2.10D . The first step is to substitute Eq. ( [Al] ) 
into ( |A3| ) and to introduce the adjoint operators ^Dj^ 
such that this can be written in the equivalent form 



X / drr-^r\T^Ur)sDX"°ir). (A4) 

/2M 



Equation ( A4 ) can be obtained from ( [A^ ) by performing 
a number of integration by parts, and the adjoint oper- 
ators as thus determined. They are found to be given 
by 



r'fCr^f-^Cr 



-3 



s = 
s = -1 

s = -2 



(A5) 



where C = d/dr* — iut. It is useful to note that the 
operator adj oint to £ is £t — —r'^fC{r'^f)^^. With this 



known, Eq. (A5) can be obtained directly from (A2). 

The next step is to relate the functions R^f°{r) to 
the Regge- Wheeler functions X^f"{r), introduced in 
Eq. ( |2.7| ). A straightforward calculation (see Ref. |3^] 
for details) shows that these are related by the Chan- 
drasekhar transformation |47 



X, 



H,oo 



(A6) 



where 



H ^ 16(l-2»Mc.)(l-4»Af^) 



-\/' — 



(A7) 



and where 



a = u;r''fCf-'Cr 



(A8) 



is a second-order differential operator. Equation (A6) 
implies 



= -4(1 - 2iMLu){l - 4iMLu){MLjyA'^f, (A9) 

the constant A^^^ was introduced in Eq. (E 

In order to re-express Z^^^ as in Eq. ( ^.1C| ) we must 
operate with sDl^C^ on the Regge- Wheeler functions. 
These operations involve successive applications of £, 
and the Regge- Wheeler equation ( |2.6| ) is substituted as 
often as necessary to simplify the expressions. To carry 
out these manipulations efficiently, it is useful to first de- 
velop the algebra of the C and C operators. We now list 
the most useful rules. 

When acting on any solution X^i{r) of the Regge- 
Whccler equation, two successive applications of £ or £ 
reduce to: 



CC = 2iijjC + V{r), 
CC^CC^ V{r), 
CC^ -2iujC + V{r) 



(AlO) 



Here , V{ r] is the Regge- Wheeler potential, introduced in 
Eq. ( |2.6| ). When carrying out the calculations it is also 
useful to invoke the commutation relations [£,g(r)] = 
[£,g(r)] — fdg/dr, where g is any function of r. 
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Using these rules the calculat ions a re straightforward. 
The final answer takes the form ( ^.10| ) if we define sTuie = 
co-'r-^DlC^, and = -&x"r Ki- 
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